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Sudden cardiac death (SCD), which is primarily caused by 
lethal heart disorders resulting in structural and arrhyth- 
mogenic abnormalities, is one of the prevalent modes of 
death in most developed countries. Myocardial ischemia, 
mainly due to coronary artery disease, is the most com- 
mon type of heart disease leading to SCD. However, post- 
mortem diagnosis of SCD is frequently complicated by 
obscure histological evidence. Here, we show that certain 
mRNA species, namely those encoding hemoglobin A1/2 
and B (Hba1/2 and Hbb, respectively) as well as pyruvate 
dehydrogenase kinase 4 (Pdk4), exhibit distinct postmor- 
tem expression patterns in the left ventricular free wall of 
SCD subjects when compared with their expression pat- 
terns in the corresponding tissues from control subjects 
with non-cardiac causes of death. Hba1/2 and Hbb mRNA 
expression levels were higher in ischemic SCD cases with 
acute myocardial infarction or ischemic heart disease 
without recent infarction, and even in cardiac death sub- 
jects without apparent pathological signs of heart injuries, 
than control subjects. By contrast, Pdk4 mRNA was ex- 
pressed at lower levels in SCD subjects. In conclusion, we 
found that altered myocardial Hba1/2, Hbb, and Pdk4 
mRNA expression patterns can be employed as molecular 
signatures of fatal cardiac dysfunction to forensically im- 
plicate SCD as the primary cause of death. 



Division of Forensic Medicine, National Forensic Service, Seoul 1 58-707, 
Korea, ' Department of Legal Medicine, ^Department of Anatomy and 
Institute of Human Genetics, College of Medicine, Korea University, 
Seoul 136-705, Korea, ^These authors contributed equally to this worl<. 
'Correspondence: sjsme@korea.kr 

Received 19 November, 2013; revised 15 January, 2014; accepted 17 
January, 2014; published online 19 March, 2014 

Keywords: heart failure, hemoglobin, postmortem RNA expression, 
pyruvate dehydrogenase l<inase, sudden cardiac death 



INTRODUCTION 

Sudden cardiac death (SCD) is a prevalent mode of death in 
most developed countries, witli an estimated overall incidence 
of 1 to 2 per 1 000 population per year. Over the last decade, 
significant advances have been made in understanding the 
clinical and genetic basis of SCD, which is primarily caused by 
heart disorders with structural and arrhythmogenic abnormali- 
ties. In forensic cases, ischemic heart disease, caused mainly 
by coronary artery disease, is the leading cause of SCD, ac- 
counting for more than half of all cases. Both acute myocardial 
infarction and ischemic coronary disease without recent infarc- 
tion are major causes of ischemic SCD. Dilated and hyper- 
trophic cardiomyopathies as well as congenital heart diseases 
such as congenital long QT syndrome, catecholaminergic po- 
lymorphic ventricular tachycardia, and Brugada syndrome can 
also result in fatal arrhythmia and cardiac death. However, 
forensic diagnosis of cardiac dysfunction associated with the 
sudden onset of arrhythmia or heart failure in autopsy samples 
is frequently complicated because it is difficult to interpret the 
severity of or determine the difference between preexisting and 
terminal dysfunctions based on conventional morphological 
findings (Oliva et al., 201 1 ). 

Postmortem RNA expression profiling, which is an emerging 
field of forensic medical science, may provide alternative tools 
for forensic diagnosis of SCD in autopsy cases. Indeed, exten- 
sive studies are currently been carried out to address key is- 
sues in forensic science regarding the cause and process of 
death, postmortem interval, and pathophysiological conditions 
of diseases and injuries (Maeda et al., 2010; Vennemann and 
Koppelkamm, 2010a). Distinct RNA compositions in autopsy 
materials have been proposed to imply certain causes and 
circumstances of death (Chen et al., 2012; Chung et al., 2012; 
Ikematsu et al., 2005; Miyazato et al., 2012; Takahashi et al., 
2009; Zhao et al., 2006; 2009; Zhu et al., 2008). Notably, Mae- 
da and colleagues proposed that left ventricular (LV) expres- 
sion profiles of several mRNA species encoding hypoxia-indu- 
cible factor 1a and natriuretic peptides were relatively higher in 
autopsy materials of cardiac death subjects, particularly those 
from subjects with ischemic heart injuries (Chen et al., 2012; 
Zhu et al., 2008). However, the application of molecular pathol- 
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Table 1 . Subject information 



Group 




Age (years) 


BMI 


Note 


Size/Sex — 


Range (mean ± SEM) 


(mean + SEM) 


TD 


n = 16/M 


25-73 (45.47 ± 2.84) 


23.69 + 0.85 


No injuries in cardiac tissues 


IH 


n = 19/M 


40-69 (52.68 ± 1 .44) 


26.91 + 2.49 


Old infarction in some cases (10/19) 


Ml 


n = 15/M 


40-63 (54.60 ± 2.00) 


25.22 + 0.53 


Calcification/thrombosis (15/15) 


Nl 


n = 20/M 


22-72 (47.10 ±3.23) 


23.27 + 0.86 


No ischemic pathology 



TD, Traumatic death; IH, Ischemic heart disease without recent infarction; Ml, Acute myocardial infarction; Nl, Non-ischemic sudden cardiac death 



ogy to diagnose unexpected cardiac death is still challenging 
because of the lack of well-established forensic RNA biomark- 
ers of fatal cardiac dysfunction. 

Our goal in this study was to identify novel molecular signa- 
tures of cardiac dysfunction after death for the forensic diagno- 
sis of SCD. For this purpose, we explored the pattern of ex- 
pression of differentially expressed genes (DEGs) associated 
with both ischemic and non-ischemic SCD and compared it 
with patterns found in non-cardiac death control subjects by 
examining genome-wide postmortem RNA expression in au- 
topsy materials from the LV free wall. 

MATERIALS AND METHODS 

Subjects and tissue preparation 

Non-infarcted myocardial tissue from the LV free wall was col- 
lected from selected medicolegal autopsy cases within 5 days 
postmortem with the help of the National Forensic Service of 
Korea. The causes of death examined in the present study 
were traumatic death without injuries to the heart (TD), ischem- 
ic heart disease (IH) with or without old infarction, acute myo- 
cardial infarction (Ml), and non-ischemic SCD (Nl), as classified 
by routine morphological, biochemical, and toxicological find- 
ings (Table 1). Dissected tissues were treated immediately with 
RNAIaterT" solution (QIAGEN GmbH, Germany) according to 
the manufacturer's instructions and stored at -70°C until use. 

RNA isolation and microarray analyses 

Total RNA was isolated from 20-50 mg tissue using the RNea- 
sy mini kit (QIAGEN) and then quantified by spectrophotometry. 
Gene expression profiling by microarray analysis was per- 
formed to explore SCD-associated human RNA transcripts as 
described previously with minor modifications (Chung et al., 
2012). Group-pooled cardiac RNA samples representing TD, 
CC, Ml, and Nl cases were subjected to microarray analysis. 
Two hundred nanogram aliquots of total RNA were amplified 
and labeled. Five micrograms of labeled cDNA was then hybri- 
dized to GeneChip Human Gene 1 .0 ST arrays (Affymetrix Inc., 
USA), and the arrays were then washed, stained, and scanned 
according to the protocol described in the Affymetrix GeneChip 
Expression Analysis Manual (Affymetrix Inc.). Signal intensities 
of each gene of interest were then normalized to those of glyce- 
raldehyde 3-phosphate dehydrogenase (Gapdh). 

Real-time quantitative reverse transcription-polymerase 
chain reaction (qRT-PCR) 

qRT-PCR reactions were performed as described previously 
with minor modifications (Chung et al., 2012). Relative mRNA 
levels were deduced according to a previous report (PfaffI et al., 
2002). Briefly, 1 ng of each total RNA sample was reverse-trans- 
cribed using MMLV reverse transcriptase (Promega, USA) by 
the random priming method. Then, aliquots of cDNA were sub- 



jected to quantitative real-time PCR in the presence of SYBR 
Green I (Life Technologies Corp., USA). Reference cDNA was 
prepared by pooling aliquots of each sample and its three-fold 
serial dilutions were utilized to calibrate PCR reactions as well 
as to construct a regression curve for calculating PCR efficien- 
cies of each cycle. Hemoglobin a1/a2 {Hba1/2), hemoglobin p 
{Hbb), and pyruvate dehydrogenase kinase 4 (Pdk4) mRNA 
expression levels were normalized to those of Gapdh. Because 
Hba1 and Hba2 mRNA share identical coding sequences, 
these two isoforms were not distinguished in the present study. 
All PCR products were cloned into the pGEM-T vector (Prome- 
ga), and their sequence identities were confirmed by chain ter- 
mination sequencing. Primer sequences used for real-time 
qRT-PCR were as follows: Hba1/2 up, 5'- TTA AGG GCC ACG 
GCA AGA AG-3'; Hba1/2 dn, 5'- GCA GTG GCT TAG GAG 
CTT GA-3'; Hbb up, 5'- GTG AAC GTG GAT GAA GTT GG-3'; 
Hbb dn, 5'- ACT TTC TTG CCA TGA GCC TT-3'; Pdk4 up, 5'- 
GCT GAT GAA CCA GCA CAT TC-3'; Pdk4 dn, 5'- CAC TCA 
AAG GCA TCT TGG AC-3'; Gapdh up, 5'- ATG TTC GTC ATG 
GGT GTG AA-3'; Gapdh dn, 5'- GGT GCT AAG CAG TTG 
GTG GT-3'. PCR efficiencies (E) for each set of primers from 
three independent measurements were as follows (in mean ± 
SEM): ^Hbai/2 = 2.0526 ± 0.001 2; Ehi* = 2.0349 ± 0.01 01 ; Ep*4 
= 2.0827 ± 0.0053; Eaapan = 2.0557 ± 0.0069 

Statistical analysis 

Data were statistically evaluated by Student's f-test or one-way 
analysis of variance (ANOVA) followed by post-hoc comparison 
by the Newman-Keuls test. Statistical significance was set at p 

< 0.05. 

RESULTS 

Cardiac and coronary pathology 

We first examined pathological features of the subjects used in 
the present study to justify our medicolegal classification of the 
cases. In comparison with TD as controls, SCD cases showed 
a significant increment in heart weights (Fig. 1A; F(3,6i)= 9.21, p 

< 0.01 by one-way ANOVA; p < 0.01 for IH and Ml, and p < 
0.05 for Nl vs. TD by post-hoc comparison) as well as LV wall 
thickening (Fig. 1B; F(3,5i)= 5.39, p < 0.01 by ANOVA; p < 0.01 
for IH, and p < 0.05 for Ml and Nl vs. TD). SCD subjects with 
ischemic heart injuries (IH and Ml) exhibited more severe car- 
diac hypertrophy features than Nl cases. Three SCD groups 
were clearly distinguishable based on the extent of main occlu- 
sions in the left anterior descending artery (LAD; Fig. 1C) and 
right coronary artery (RCA; Fig. ID) (for LAD, F(3,63)= 61.07, p 

< 0.01 by ANOVA and p < 0.01 for IH/MI vs. TD/NI by post-hoc 
comparison; for RCA, Fp, 63) = 27.66, p < 0.01 by ANOVA and p 

< 0.01 for IH/MI vs. TD/NI). Moderate coronary stenosis was 
found in IH subjects (60.26 ± 7.14% for LAD; 52.1 1 ± 5.55% for 
RCA), while severe or complete blockade was observed in the 
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Table 2. Differentially expressed genes 



Gene symbol 


Gene accession — 


IH 


Fold cfiange 
Ml 


Nl 


1 AIJI COOIUI 1 ICVCIO 

{% Of GAPDH) 


Hba1/2* 


NM_000558/NM_00051 7 


1.84 


2.31 


2.19 


57.47 


Hbb 


NM_000518 


1.86 


2.25 


2.24 


75.53 


Areg 


NM_001657 


-2.29 


-2.48 


-1.99 


1.33 


B3gnt5 


NM_032047 


-2.94 


-1.68 


-2.09 


1.72 


Fcn3 


NM_003665 


-2.01 


-2.12 


-1.83 


3.32 


Nr4a2 


NM_006186 


-3.01 


-1.64 


-2.72 


3.51 


Pdk4 


NM_002612 


-2.92 


-3.06 


-6.80 


21.64 



"Hbal and Hba2 are considered to be as single gene here because they have identical coding sequences and can therefore not be distinguished 
based on microarray or RT-PCR analyses. 
"Average expression levels of the four groups for the indicated genes relative to that of Gapdh. 




TD IH Ml Nl TD IH Ml Nl 



C D 




TD IH Ml Nl TD IH Ml Nl 



Fig. 1. Cardiac and coronary features according to different causes 
of death. Heart weight (A) and the thickness of the left ventricular 
free-wall (B) were measured in subjects who died due to traumatic 
death (TD), ischemic heart disease (IH), acute myocardial infarction 
(Ml), or non-ischemic sudden cardiac death (Nl), and are expressed 
as means ± SEM. Extent of occlusion at the site where a major 
stenosis was found in the left anterior descending arteries (LCA; C) 
or the right coronary arteries (RCA; D) is expressed as the mean ± 
SEM (%). Data were statistically evaluated by one-way ANOVA 
followed by the Newman-Keuls test (*p < 0.05 and ** p < 0.01 vs. 
TD). 



Ml cases (93.33 ± 1 .52% for LAD, p < 0.01 for IH vs. Ml; 81 .67 
± 7.77% for RCA, p < 0.01 for IH vs. Ml). In contrast to SCD 
subjects with ischemic heart injuries, the extent of LAD and 
RCA occlusion in Nl subjects was essentially the same as that 
observed in controls (usually less than 20%). 

Identification of SCD-associated differentially expressed 
genes (DEGs) 

To identify mRNA species associated with SCD, we carried out 
microarray analysis on postmortem RNA samples derived from 



LV myocardial tissues. In our initial approach, we explored 
individual genes differentially expressed between controls (TD) 
and SCD subjects categorized into the three subgroups of IH, 
Ml, and Nl. Out of 13,288 RNA species that showed differential 
expression between control and SCD tissues, 85 (IH), 517 (Ml), 
and 514 (Nl) DEGs showed a more than 1.5-fold difference in 
expression compared with the TD controls. Among these, eight 
gene transcripts were common to all subgroups of SCD (Table 
2). mRNA expression levels of hemoglobin a1/a2 (Hba1/2) and 
hemoglobin p (Hbb), which are components of the adult he- 
moglobin tetramer, were higher in all SCD cases than controls, 
whereas several gene transcripts such as pyruvate dehydroge- 
nase kinase 4 (Pdk4), nuclear receptor 4, subgroup A, member 
2 {Nr4a2), and ficolin 3 {Fcn3} showed lower expression in 
SCD cases. Considering relative expression levels (expressed 
as % of Gapdh as shown in Table 2) as well as fold-differences 
between SCD and control cases, we focused mainly on Hba1/2, 
Hbb, and Pdl<4 mRNAs to examine whether alterations in levels 
of these mRNAs in individual autopsy materials could suggest 
fatal cardiac dysfunction leading to SCD. 



Hba1/2, Hbb, and Pdk4 mRNA profiles in individual 
subjects with known causes of death 

In the next set of experiments, we compared cardiac Hba1/2, 
Hbb, and Pdk4 mRNA profiles in individual subjects with three 
different modes of SCD to those in TD cases using qRT-PCR 
results. Hba1/2 and Hbb mRNA levels tended to be higher in all 
SCD-related groups than in the TD group by approximately 4- 
fold (Fig. 2A; for Hba1/2: F(3,ee)= 2.33, p = 0.082 and for Hbb: 
F(3, 66) = 2.35, p = 0.081 by ANOVA; p < 0.05 vs. TD by Stu- 
dent's f-test for each SCD-related group). In contrast, Pd/c4 
mRNA levels were significantly lower in the SCD-related sam- 
ples (F(3,66)= 4.60, p < 0.01 by ANOVA; p < 0.05 for TD vs. Nl 
and p < 0.01 vs. TD for IH or Ml cases by Neuman-Keuls post- 
hoc comparison). Although these mRNA species exhibited 
group differences according to cause of death, large variations 
within groups were found. Box-and-whisker diagrams for diffe- 
rentially expressed mRNA transcripts show that more than one- 
half of measured values in a group overlapped with values in 
other groups, indicating poor discrimination (Fig. 2B). In con- 
trast to myocardial tissues, Hba1/2, Hbb, and Pdk4 mRNA levels 
in the occipital cortices of the brain were not different between 
TD and SCD subjects (Fig. 3). Brain tissues are known to be 
highly vulnerable to hypoxic/anoxic states leading to cellular 
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Fig. 2. Hba1/2, Hbb, and Pdl<4 nnRNA expression in the left ventri- 
cular free wall in association with different causes of death. Hba1/2, 
Hbb, and Pdk4 mRNA expression profiles were quantitatively ana- 
lyzed by real-time qPT-PCR in subjects who died due to traumatic 
death (TD), ischemic heart disease (IH), acute myocardial infarction 
(Ml), or non-ischemic sudden cardiac death (Nl). (A) Data were 
normalized to Gapdh mRNA levels and expressed as the mean ± 
SEM in arbitrary units (A.U.), where the mean values of the TD 
group were set to 1 . Data were statistically evaluated by one-way 
ANOVA followed by the Newman-Keuls test (*p < 0.05 and **p < 
0.01 vs. TD). (B) Relative mRNA levels were expressed by box- 
and-whisker diagrams with log scales. Median values are indicated 
by the lines in the middle of the boxes. 



stress responses and metabolic disturbances during the agonal 
stage as well as the acute postmortem period (Franz et al., 
2005). Therefore, our results strongly suggest that SCD-related 
alterations in Hba1/2, Hbb, and Pdk4 mRNA expression reflect 
myocardial dysfunction rather than cellular damage during 
agonal and/or postmortem periods. 

Ratios of Hba1/2 or Hbb mRNA levels relative to those of 
Pdk4 as potential forensic biomarl<ers 

Despite the group differences in Hba1/2, Hbb, and Pdl<4 mRNA 
expression depending on cause of death, large variations within 
the group were found in all selected genes as shown in Fig. 2B. 
These features suggest that Hba1/2, Hbb, or Pdk4 mRNA ex- 
pression relative to that of a housekeeping gene, such as Gapdh 
employed in the present study, may not singly provide practical 
information by itself. However, concurrent alterations in the 
relative composition of the marker mRNA species could be 
used as a valuable index to better discriminate among causes 
of death, as we demonstrated in our previous study of mRNA 
transcripts encoding heat shock proteins (Chung et al., 2012). 
To test this idea, we bi-directionally plotted relative expression 
levels of Hba1/2 or Hbb mRNA normalized by Gapdh in relation 
to those of Pdk4. All SCD-related groups occupied positions 
distinct from those of TD subjects in the plot (Fig. 4A). We 
therefore compared ratios of Hba1/2 or Hbb mRNA abundance 
to those of Pdk4. As shown in Fig. 4B, the relative ratio of 
Hba1/2- or Hbb-to-Pdk4 expression levels was higher in all SCD- 
related groups, including Nl subjects, compared to TD controls 
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Fig. 3. Hba1/2, Hbb, and Pdk4 mRNA expression in the occipital 
lobes of cerebral cortices. Hba1/2, Hbb, and Pdk4 mRNA expres- 
sion profiles were analyzed quantitatively by real-time qRT-PCR in 
subjects who died due to traumatic death (TD) or sudden cardiac 
death (SCD) including both ischemic and non-ischemic heart inju- 
ries. (A) Data were normalized to Gapdh mRNA levels and ex- 
pressed as the mean ± SEM in arbitrary units (A.U.), where the 
mean values of the TD group were set to 1. (B) Relative mRNA 
levels are shown in box-and-whisker diagrams with log scales. 
Median values are indicated by the lines in the middle of the boxes. 



(for the Hba1/2-to-Pdk4 ratio, Fp, 66)= 4.3180, p < 0.01 by one- 
way ANOVA; for Hbb-to-Pdk4 ratio, F(3,66)= 2.8520, p < 0.05). 
In box-and-whisker diagrams of these ratios, overlapping sec- 
tions were reduced by approximately -25%, indicating better 
resolution between TD and SCD subjects (Fig. 4C). Taken to- 
gether, our findings suggest that concurrent alterations in Hba1/2, 
Hbb, and Pdk4 expression may be good molecular indices for 
postmortem diagnosis of sudden unexpected cardiac death 
due to a wide range of causes. 

DISCUSSION 

In the present study, we aimed to identify LV myocardium- 
enriched mRNA species detectable in postmortem heart tis- 
sues with expression profiles associated with fatal cardiac dys- 
functions leading to SCD. For this purpose, we compared ge- 
nome-wide mRNA expression profiles in LV free wall tissues 
from subjects with three different modes of SCD and compared 
these profiles with those obtained for TD cases as a control 
group. We found that patterns of Hba1/2, Hbb, and Pdk4 
mRNA expression in SCD cases were distinct from that ob- 
served in TD subjects, but common among subgroups of SCD 
subjects. 

It should be noted that the concept of molecular or genetic 
autopsy was first established in postmortem investigation of 
SCD. Early molecular studies of SCD cases included detection 
of viral genomes in inflammatory cardiomyopathies as well as 
gene mutational analysis of several hereditary heart diseases 
(Basso et al., 2001; Tester and Ackerman 2006). Along with 
recent advances in molecular genetics, several genetic defects 
in cardiac ion channels such as KCNQ1, KCNH2, KCNE1/2, 
and SCN5A have been shown to be associated with certain 
cases of heritable and acquired heart malfunctions due to ven- 
tricular arrhythmias and sudden death. 

However, it is also believed that RNA expression analysis of 
autopsy tissues from pathologic and non-pathologic hearts 
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Fig. 4. Pairwise analyses of relative mRNA abun- 
dances of Hba or Hbb in relation to Pdk4. (A) 
Pairwise analyses of Hba1/2 or Hbb mRNA levels 
in comparison with those of Pdk4. Normalized 
mRNA levels of Hba1/2 (upper) or Hbb (lower) 
relative to Gapdh were bi-directionally plotted in 
relation to those of Pdk4 on a logarithmic scale. 
Each symbol and error bar indicates the mean + 
SEM. (B) Ratio of relative Hba1/2 (upper) or Hbb 
(lower) mRNA level to that of Pdk4 in each case 
were calculated and expressed as means + SEM. 
Data were statistically evaluated by one-way 
ANOVA followed by the Newman-Keuls test (*p < 
0.05 and **p < 0.01 vs. TD). (C) The log ratios of 
Hba1/2 or Hbb to Pdk4 were expressed as box- 
and-whisker diagrams. Median values are indi- 
cated by the lines in the middle of the boxes. 



might be useful for postmortem investigation of SCD cases. 
Both qualitative and quantitative features of cardiac mRNA 
expression have been examined. For example, an increase in 
non-functional splice variants of the ScnSa gene transcript was 
show/n to be associated with heart failure due to impairment of 
the cardiac Na* current along with a reduction in NaV1 .5 chan- 
nel protein levels (Shang et al., 2007). In addition, Maeda and 
colleagues showed that several hypoxia-responsive gene tran- 
scripts encoding hypoxia-inducible factor-1a, erythropoietin, 
and vascular endothelial growth factor were locally increased in 
cardiomyocytes from acute myocardial infarction cases com- 
pared with asphyxiation cases, but not in recurrent myocardial 
infarction and uncharacterized SCD cases (Zhu et al., 2008). In 
a more recent study, these authors also reported that mRNA 
expression of atrial and brain natriuretic peptides in the LV free 
wall tended to be higher in various modes of cardiac death 
including acute ischemic heart disease, chronic congestive 
heart disease, and right ventricular cardiomyopathy than in LV 
tissues from acute fatal bleeding and/or mechanical asphyxia- 
tion cases (Chen et al., 2012). Therefore, it is crucial as well as 
practically useful to explore sets of myocardial RNA transcripts 
closely and commonly related to various modes of heart failure 
for postmortem diagnoses of SCD cases, even those without 
apparent anatomical signatures. 

Hba1/2 and Hbb mRNAs have been believed to be predomi- 
nantly expressed in cells of erythroid lineage. Recent studies, 
however, reported HBA1/2 and HBB expression in non-ery- 
throcyte cells such as neurons (Richter et al., 2009), alveolar 
epithelial cells (Newton et al., 2006), endometrial cells (Dassen 
et al., 2008), hepatocytes (Liu et al., 2011), as well as cardiac 
cells (Arab et al., 2007; Kittleson et al., 2005; Solskov et al., 
2012). In particular. Hare and colleagues showed that Hba1/2 
and Hbb mRNA levels were upregulated in myocardial tissues 
biopsied from both ischemic and non-ischemic cardiomyopathy 
patients (Kittleson et al., 2005). Although the biological function 
of cellular hemoglobin in non-erythroid cells requires elucidation, 
cellular expression of these hemoglobin subunits appears to be 
induced in response to increased levels of intracellular reactive 
oxygen species (ROS), and they may play a protective role 
against oxidative stress (Arab et al., 2007; Li et al., 2013; Liu et 
al., 2011). In this context, it is noteworthy that oxidative stress, 
characterized as increased generation of mitochondrial ROS 
production in cardiac cells, is associated with various modes of 



heart failure (Belch et al., 1991; Hill and Singal, 1997; Matsu- 
shima et al., 2006). Prolonged ROS production can result in a 
variety of cardiac dysfunctions and damage the failing heart; it 
can impair cardiac contractility, activate a variety of hypertroph- 
ic signaling cascades and transcriptional networks, and directly 
mediate apoptosis (Tsutsui et al., 201 1). It is therefore plausible 
that Hba1/2 and Hbb mRNA expression might be induced by 
oxidative stress during myocardial dysfunction leading to car- 
diac death. 

PDK4, a mitochondrial protein, is known to play an important 
regulatory role in maintaining metabolic flexibility in cardiomyo- 
cytes by causing a switch from predominant oxidation of fatty 
acids to oxidation of glucose by phosphorylating pyruvate de- 
hydrogenase (PDH). Loss of metabolic flexibility, such as that in 
diabetes mellitus, is closely related to the onset of cardiovascu- 
lar diseases and heart failure (Alvarez-Guardia et al., 2010; 
Taegtmeyer et al., 2002; Zhao et al., 2008). It was recently found 
that a splicing mutant of the Pdk4 gene lacking exon 1 1 , which 
is a prerequisite for PDK4 protein activity, underlay the patho- 
physiological phenotypes in a Doberman pinscher model of 
human dilated cardiomyopathy (Meurs et al., 2012; Wynn et al., 
2008). More importantly, a systemic analysis to discover gene 
networks in the LV myocardium revealed that peroxisome proli- 
ferator-activated receptor gamma coactivator-la (PGC-1a) tar- 
get gene down-regulation may serve as a key signature of a 
failing human heart (Sihag et al., 2009). Impaired PGC-1a ac- 
tivity was shown to lead directly to a reduction in PDK4 expres- 
sion, resulting in an increment in glucose oxidation (Alvarez- 
Guardia et al., 2010). Therefore, it is tempting to speculate that 
down-regulation of cardiac Pdk4 mRNA expression, either for 
genetic or environmental reasons, might be related to suscepti- 
bility to cardiac dysfunction. This possibility needs to be further 
examined. 

From a practical viewpoint, routine normalization methods 
using a housekeeping gene in quantitative analyses of post- 
mortem RNA expression often mislead and may obscure dif- 
ferences between cases and controls due to a variety of ante- 
mortem and postmortem factors differentially influencing integri- 
ty of each RNA species (Gonzalez-Herrera et al., 2013; Preece 
and Cairns, 2003; Zhang et al., 2013). Rather, pain/vise analys- 
es of a set of gene transcripts with concurrently changed, but 
negatively correlated under certain pathological condition, ap- 
peared to be more practical for forensic applications as we 
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previously demonstrated (Chung et al., 2012). As shown in Fig. 
4, Hba1/2 or Hbb mRNA expression was negatively correlated 
with Pdk4 mRNA expression, and paired analyses of two 
mRNA species (Hba/b and Pdk4) provided better indices imply- 
ing cardiac death. 

In conclusion, we suggest that concurrent alterations in the 
expression of heart-enriched Hba1/2, Hbb, and Pdk4 mRNAs 
in the postmortem LV free wall may reflect a wide spectrum of 
fatal cardiac dysfunctions. Expression profiles of a few mRNA 
transcripts may not be sufficient for postmortem diagnosis of 
SCD cases, especially non-ischemic cases; further evaluation 
of heart failure-related sets/pairs of gene transcripts would 
therefore be useful to enable practical application of such DEGs 
to forensically diagnose SCD. 
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